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1. Introduction. Up to date, the analysis of air cascade family has 
been performed mainly by using a full Monte Carlo simulation, because 
of the large fluctuation of the air cascade. Recently it has become 
clear, however, that it is difficult to draw a definite conclusion about 
the interaction mechanism by using only this kind of simulation, due to 
too many assumptions in relation to the amount of information. 

On the other hand, some attempts to reproduce the original y-ra y at 
the interaction point, for example ''decascading 1 ", have also been made. 
This kind of method makes it possible to observe the interaction di- 
rectly and to analyze the data from various angles. All of these me- 
thods, however, assume the constant ER in cascade shower, where E is 
energy and R is the distance from the center of cascade shower. This 
character is all-inclusive and changes according to the energy and the 
interaction height of original y-ray. It is impossible, therefore, to 
reproduce the exact interaction height and energy by these methods. 

This work adopts a quite different way from the ER constant method; that 
is, a relative method in separating one cascade shower from others. This 
new method makes it possible to estimate the interaction height and 
energy by using information about the lateral spread of the cascade 
shower . 

2. New method. Fig.l (a) shows the correlation between E and R of y- 
rays in a simulated family event. Open circle and closed circle repre- 
sent the Y-ray produced by one original Y-ray and other original y-rays, 
respectively. In this figure, we can find an energy gap at the boun- 
dary between groups of cascade shower Y-rays. The new method is noth- 
ing but a way to classify one cascade shower group from others by using 
of this energy gap. 

We make a classification of each cascade shower according to the fol- 
lowing procedure. 


* 





Fig.l (a) Correlation between E and R in a simulated family event. 
Open and closed circle represent the y-ray produced by one original y- 
ray and othe y-rays, respectively, (b) Transformed figure of (a) in 
order to represent the energy gap symmetrically . 
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1) Take the i-th y- r ay as a temporary center of a cascade shower group. 
The selection of the i-th y-ray is performed according to size of ener- 
gy, because a high energy y-ray should be located near the true center. 

2) Calculate the following quantity Sij representing a spread between 
the i-th y-ray of a temporary center and all other Y-rays ( j* l-n ) , 

Sij = (Ei+5)Rijv^j/l5 (l) 

where Rij is a length between the i-th and j-th y-rays. As is seen in 
Fig.l (a), a straight line of Rv^E = constant has a mean inclination be- 
tween the line of ER=constant and R=constant. 

3) Caluculate the energy density for the following quantity W j , 

Wj = 1-exp (-Sj ) (2) 

for fixed i. Wj ranges from 0 to 1 and is effective to search for an 

energy gap impartially among various cases. The correlation between E 
and W is shown in Fig. 1 (b) for the same sample as in (a). We can see 
a symmetrical valley of energy in this Figure. 

The energy density is defined as 

Pj = ^ Ej*Ej + i*Ej +2 /(Wj +2 -Wj) (3) 

The boundary is set up the position Wmin where Pj has the minimum value. 
b) Select candidates, which satisfy the following criteria, as members 
of the same cascade shower with the i-th y-ray. 

Wj < Wmin and Ej < Ei 

5) Proceed to the next y-ray to be set as a temporary center. The y- 
rays which have been already selected in (^) are excluded as new centers. 

6) When a y-ray is selected as a candidate of more than two groups, it 
is decided to belong the group which gives the smallest distance from 
its center. 

3. Estimation of energy and production height. The energy and produc- 
tion height of y-ray are estimated by using the mean spread of cluster- 
ized y-ray as members of a cascade shower. The mean lateral spread 
<R> is defined as, 

<R> = ZR/N (cm) (k) 

where R is a distance from the energy weighted center and N is the num- 
ber of y-rays in a cluster. Close investigation of the simulation data 


drew the estimation formula of energy Ees 
Mt. kanbara height (5500m) as follows. 

and 

interaction height Hes at 


Ees = EE •10 2,5<R> 

(TeV) 

( 5 ) 


He s == 17000 *<R> 1,5 

(m) 

(6) 

w 


b. Check of new analysis method. Reproducibility by the new method 
was checked by using artificial data in a full Monte Carlo simulation. 
For the purpose of checking, only the standard model simulation was nec- 
essary. An outline of the simulation 2 is as follows. 

(1) Proton primary with energy spectrum as 

I ( E ) dE « E" 6-1 dE ( ) 

(2) Particle production spectrum is scaling law 
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* 


f (x)dX a (l-X) iS /XdX 

where X=E/Eo, E and Eo are secondary and primary energy, respec- 

tively. 

(3) Interaction mean free path 

X = / 80 g/cm 2 for proton 

\ 96 g/cm 2 for secondary pion 

( b ) Mean transverse momentum 

= / 330 MeV/c for secondary particle 

\ 500 MeV/c for leading particle 

For the investigation of large Pt, the double value <Pt>=660 MeV/c was 
also used. 

All of the clusterized groups have not always one to one corres- 
pondence to the expected one. In the case of no one to one correspond- 
ence, the original Y-ray which gives the highest energy ratio within the 
cluster is selected as corresponding y-ray and used for checking. 

Fig. 2 (a) shows the distribution of Nes/No, where Nes is the number 
of y-rays in an actual clusterized group and No is the expected number 
in the corresponding cascade shower. The peak at Nes/N 0 = 2 is due to 
the case of Nq = 1. Figure shows fairly strong agreement and the vali- 
dity of the new method was confirmed. 

Fig. 2 (b) and (c) show the comparison of estimated energy Ees and 
estimated interaction height Hes of original y-ray with expected true 
ones Eo and Ho, respectively. Fig. (c) for interaction height shows 
broad distribution. This is due to mainly the fluctuation in the tra- 
veling distance of the original y-ray from the interaction point to the 
first pair creation point. We need the statistical analysis for the 
interaction height. 



Fig. 2. Comparison of the estimated value by new method with the 
expected one (a) Number of y-rays in a clusterized shower (b) 
energy and (c) interaction height of original Y~ ra y 


Fig. 3 shows the reproducibilty of the transverse momentum Pt distri- 
bution, (a) for normal Pt and (b) for double Pt. The full and broken 
lines represent the distribusions of the estimated Pt and of the expect- 
ed Pt of all Y-rays at the interaction point. The estimated and expec- 
ted distribtion does not necessarily coincide because of the detection 
bias and the error of estimated value. The positions of peak, however, 
coincide with each other, and the estimated distribution for normal Pt 
and for large Pt reproduce just a double difference. It is worth notice 
that the estimated Pt distribution has a large error but is independent 
of the interaction model. 
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Fig . h The interaction 
height distribution, 
compared with the simu- 
lated data(broken line) 



Fig. 5 The lateral 
spread distribution of 
original y-ra y. Smooth- 
ed curves represent the 
simulated data for 
normal. Pt and large Pt. 



Fig . 6 The transverse 
momentum Pt distribu- 
tion of reproduced ori- 
ginal y~ray, compared 
with the simulated dis- 
tribution for normal Pt 
(<Pt>=330 MeV/c) and 
large Pt (660 MeV/c) 


Fig. 3 The transverse momentum Pt distri- 
bution of reproduced y-ray for (a) normal 
Pt (b) double Pt, compared with the dist- 
ribution of all y-rays produced by the 
interaction (broken line). 

5. Results. In Fig. 4, is shown the result 
of the interaction height of the reproduced 
original Y-ray by the experimental data of 
the Mt. Kambara emulsion chamber experiment, 
compared with the result of simulated data 
(broken line). The experimental data does 
not contradict the prediction of the simula- 
tion . 

Fig. 5 shows the lateral spread R distri- 
bution, where R is the distance between an 
energy weighted center of a family and of 
Y-rays in a cluster. Smoothed curves re- 
present the distribution of the simulated 
data in the case of normal Pt and double Pt. 

Fig. 6 shows the transverse momentum Pt 
distribution by Mt.Kanbara experiment, 
compared with the simulated data. Smoothed 
curves in the figure are also the simulated 
distribution for normal Pt and double Pt. 

Though the statistics of experimental 
data are not enough. to reach a definite con- 
clusion, the both results in Fig. 5 and Fig. 

6 prefer large Pt ( <Pt> - 600 MeV/c ) 
rather than normal Pt. 
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